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The scattering of two Dashen-Hasslacher-Neveu (DHN) baryons in the large A*' Gross-Neveu model 
is solved exactly using the relativistic, time-dependent Hartree-Fock approach. Unlike the special 
case of kink-antikink scattering, the scattering of DHN baryons is sensitive to the back-reaction 
of the fermions bound inside the baryons. Correspondingly the solution is much more complicated 
, than the kink-antikink scattering solutions, which can be expressed in terms of sinh-Gordon solitons. 

■ Nevertheless, we present a simple ansatz form that gives closed analytic expressions for both the 

space-time dependent mean fields and the Dirac spinors for all continuum and bound states. The 
solution can also be applied to the scattering of polarons and solitons in conducting polymers. 
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I. INTRODUCTION 



The simplest version of the Gross-Neveu (GN) model Q is a 1-1-1 dimensional model quantum field theory of N 
^ • species of massless, self-interacting Dirac fermions with Lagrangian 

^ : 2 

fe=l \fc=l / 

. We consider this model in the 't Hooft limit, N — > oo, with Ng^ ^ constant, where semiclassical methods become 
QQ ' exact. In spite of its apparent simplicity, the model ([1]) gives rise to quite non-trivial physical phenomena of interest 
QQ , for strong interaction particle physics as well as condensed matter physics. This has been established in works spread 
' out over several decades, in an effort which is still ongoing, as we shall see. The original work from 1974 was motivated 
, in part by the discovery of asymptotic freedom in quantum chromodynamics (QCD) a property which the GN 

model shares. The main focus in Ref. [l| was on spontaneous breaking of the discrete chiral symmetry — >■ 75^/', 
T-H \ fermion mass generation via dimensional transmutation and the scalar a meson. Soon afterwards, multifermion bound 
states (baryons) were found by Dashen, Hasslachcr and Neveu (DHN) in a study of the associated gap equation 
The gap equation solution for the static DHN baryon scalar potential is 
• ^ . 

\ <5'dhn(2;) — 1 + y tanh(j/a; — -arctanh(?/)) — y tanh(?/x + -arctanh(y)) (2) 



where the y parameter satisfies < y < 1. DHN baryons have a valence bound state which can be filled with up 
to N fermions. In the large N limit, the filling fraction v = n/N becomes a continuous parameter, so that there 
is in fact a whole one-parameter family of baryons. For self-consistency in the gap equation, the filling fraction v 
is related to y by j/ = sin(7r://2). The DHN baryons span the region from a weakly bound, non-relativistic state at 
small filling to the ultrarelativistic limit of a decoupled kink and antikink at complete filling. At large filling fraction, 
the DHN baryon looks like a bound kink-antikink molecule, and it is stable not because of topology but because 
of a balance between the kink-antikink interaction and the effect of the fermions bound to the kink and antikink. 
This leads to a direct relation between the baryon size (the distance between the kink and antikink) and the fermion 
filling fraction, giving a beautiful example of dynamical stability as well as of the Jackiw-Rebbi mechanism of fermion 
modes bound to localized defects such as kinks [H . This " anatomy" of the DHN baryon is sketched in Figure [T] A 
particularly interesting special case is the baryon with non-trivial topology, the kink, which is attributed to Callan, 
Coleman, Gross and Zee (CCGZ) (cited in Q). Later, Feinberg established the complete set of static solutions to 
the large N gap equation Q, combining inverse scattering theory with resolvent techniques. The general solution 
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consists of marginally bound multibaryon configurations whose energy does not depend on the distance between 
the constituents. A common feature of all static solutions is the fact that the self-consistent scalar potentials are 
reflectionless, generalizing the Kay- Moses potentials of the Schrodinger equation to the Dirac equation Q. 




FIG. 1: A sketch of the anatomy of a DHN baryon, here shown for parameter y — 0.999999. The scalar field 5'DHN(a;) is shown 
as the solid blue line going asymptotically to 1 (in units of the dynamically generated fermion mass in the GN model). The 
fermion density, p = i/'^i/'i for the bound valence mode is shown as a black dotted line; note that the density is localized at the 
kink and antikink of the baryon scalar field. The pseudoscalar condensate, tt{x) — ijji'ys^p, for the valence state is plotted as a 
red dashed curve, while the scalar condensate, cr{x) = ipip, for the valence state is plotted as a purple dot-dashed curve (cr has 
been multiplied by a factor of 100, for visualization purposes) . The scalar condensate is nonzero inside the DHN baryon, while 
the pseudoscalar condensate is localized at the edges, on the kink and antikink, like the density p but with a change of sign. 

The baryons of the GN model also play a key role in the properties of matter at finite density and temperature. The 
preferred state of cold matter in the GN model was shown to be a crystal, consisting of an array of kinks and antikinks 
[l(j | . The fact that the single baryon potential is reflectionless translates (at finite density) into a "finite-gap" periodic 
potential, expressed analytically in terms of elliptic functions. The phase diagram as a function of temperature and 
chemical potential displays a soliton crystal phase, a massive and a massless Fermi gas phase, and has turned out to 
be much richer than orig inally thought [11 [H. 

Semiclassical methods are not restricted to static solutions. Since we are dealing with a relativistic field theory, 
we can boost any static solution to an arbitrary Lorentz frame, turning a static Hartree-Fock (HF) solution into a 
solution of the time-dependent Hartree-Fock (TDHF) approach, 

occ 

[ilbdx + 1°S) ipa = idttpa , S = -g^ ^ ipi3i>i3- (3) 

Indeed, it was anticipated already in Witten's seminal paper on baryons in the 1/N expansion in QCD that 
baryon-baryon scattering might be solved with the help of the TDHF approach at large N. In this paper we realize 
this goal explicitly for the GN model. The TDHF approach has already been exploited to compute structure functions 
of DHN baryons and CCGZ kinks, further confirming the value of the GN model as a toy model for QCD [l^. In 
addition to boosted static solutions, some truly dynamical solutions of Eqs. ^ are also known. They are harder 
to find than static HF solutions since inverse scattering theory is not developed enough for TDHF. The efforts to 
find non-trivial, time dependent mean field solutions were also pioneered by DHN who already presented a breather 
solution, a periodically oscillating (in time) multifermion bound state DHN also pointed out the possibility to 
relate the breather to the kink-antikink scattering problem by analytic continuation. This suggestion was taken up and 
elaborated in several recent works. Kink-antikink scattering was addressed in (isj . Apart from a first glimpse of the 
scattering problem of composite, relativistic objects, the solution also gave several new insights into the mathematical 
structure of the theory. A special feature of kink dynamics is the fact that the (valence) fermions do not react back on 
the solitons that are carrying them. This decoupling made it possible to formulate kink dynamics in the language of a 
well-studied, classical soliton theory, the sinh-Gordon model [1^1, and to use the known n-soliton solution of this model 
to generalize kink-antikink scattering to the case of an arbitrary number of kink-like baryons [T7| . From the point of 
view of TDHF theory, the most striking feature of all kink-antikink solutions is the fact that the scalar density of each 
single-fermion level is proportional to the full self-consistent potential S. This kind of solution where self-consistency 
holds mode- by- mode was called type I in [l^ . This feature in turn opens the way to a geometrical interpretation of 
HF solutions of the GN model in terms of embedding 2d surfaces into 3d spaces [131 ■ Self-consistency translates into 
the condition of constant mean curvature. One can then understand the relationship to minimal surfaces in AdSa 
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and hence string worldshects in simple terms, a relation which had been found via the sinh-Gordon equation [19| but 
was otherwise completely mysterious in [l^ . 

As nice as kinks are mathematically, they form only one extreme endpoint of the DHN baryon family. To complete 
the picture, we need to understand the scattering of general DHN baryons, without the restriction to kinks and 
antikinks. The solution to the scattering problem for two arbitrary DHN baryons would allow one to probe the 
degree to which the internal bound state structure is relativistic, all the way from the non-relativistic limit to the 
ultrarelativistic one. On the other hand, by choosing the velocity of the baryons, one can cover the range from 
non-relativistic to relativistic scattering in the external kinematics as well. This situation is sketched in Figure 2. 
This general problem is not easy, as evidenced by the fact that it has never been addressed in the literature so far, to 




FIG. 2: The upper figure shows the scattering of a CCGZ-kink baryon and a CCGZ-antikink baryon, as was already discussed in 
[l5l.[l7l]. In contrast, the lower figure shows the problem treated in this current paper: the scattering of two DHN baryons. The 
solid [blue] lines show the scalar potential, while the dotted [black] lines show the density of the bound valence fermions. Note 
that each asymptotic DHN baryon can have different internal structure (i.e., a different y parameter), which affects the baryon 
size and shape as well as the internal fermion densities; this will be important for understanding the scattering processes. By 
varying the structural parameters, and also the relative velocity parameter, we can probe various interesting physical regimes, 
spanning continuously from relativistic to non-relativistic. 

the best of our knowledge. No established, systematic method for solving the TDHF equations is available. Since the 
valence fermions are now expected to react back, it is also unlikely that one can map the problem onto some known 
soliton theory, like sinh-Gordon theory for the kinks. We propose to solve this problem here by a method based upon 
an ansatz. As a matter of fact, due to the self-consistency issue, we know neither the equation nor its solution to 
start with, so that guessing the right ansatz is quite a challenge. We will explain how a simultaneous ansatz for the 
scalar TDHF potential and the (bound and continuum) spinors can nevertheless be found heuristically. The solution 
of the Dirac equation and the requirement of self-consistency are sufficient to determine the unknown parameters of 
the ansatz and to establish an exact baryon-baryon scattering solution in the large N limit of the GN model. The 
fact that this procedure is successful is undoubtedly related to the integrability of the GN model at finite N pol - [2^ , 
although we are not aware of any direct path from the integrability of the finite N GN model to the large N solution 
which we find here. Perhaps more important than the specific solution found here is the fact that the ansatz can be 
generalized in a natural way to a whole class of more complicated scattering problems, also involving multibaryon 
bound states and breathers in addition to DHN baryons. 

This paper is organized as follows. In Sec. |TT1 we illustrate our method of solving the TDHF equations via ansatz 
in a simpler context, the single DHN baryon in flight. The results are needed later on, since they enter the scattering 
problem as incoming and outgoing states. We also introduce some technicalities related to lightcone variables. Sec. lIIII 
describes the solution of the scattering problem. We first exhibit the ansatz for the scalar potential, Sec. lIII Al contin- 
uum spinors. Sec. IIII Bl and bound state spinors, Sec. IIII CI constraining the parameters by asymptotic information 
whenever possible. Sec. IIIIDI presents the results for the non-trivial parameters determined by solving the Dirac 
equation. In Sec. IIII E[ we demonstrate self-consistency of the solution and compute the fermion density. Sect. IIIIFI 
generalizes the results obtained in the ccnter-of- velocity frame to an arbitrary Lorcntz frame. In Sec. IIVI we report 
on several tests of our results carried out in limiting cases where some information is already available from other 
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sources, i.e., the static limit. Sec. IIV Al the non-relativistic limit, Sec. lIVBl and the kink-antikink limit. Sec. IIV C] 
This is followed by several numerical examples and the relationship to polarons and excitons in conducting polymers 
in Sec. IIVDI Sec. |V] contains a brief summary and an outlook. An appendix summarizes the self-consistent scalar 
potential in conventional variables, complementary to the lightcone formulation used in the main text. 



II. SINGLE BARYON IN FLIGHT 

The DHN baryon (3] is a well-understood multifermion bound state of the large N GN model. Originally derived 
in the rest frame by the inverse scattering method, it can easily be boosted to any other inertial frame. We refer 
to Ref. [3] for a detailed discussion of the baryon in flight within the Dirac-TDHF approach. In this section, we 
introduce our ansatz method by pretending that we don't know the DHN solution. This illustrates how one can find a 
self-consistent solution of the Dirac-TDHF equation ([3]) with a judiciously chosen ansatz for the scalar potential and 
the spinors. Apart from this ansatz, a crucial ingredient of our approach is the systematic use of lightcone coordinates 
and variables. These greatly simplify the final results, due to their simple Lorcntz transformation properties. We 
introduce all the necessary definitions here in a familiar context. The results for the single baryon in a general Lorcntz 
frame and in the most convenient language are obviously a prerequisite for addressing the two-baryon scattering 
problem, where they appear as incoming and outgoing asymptotic states. 



A. Solution of the Dirac-TDHF equation by ansatz 

The key quantity for any Dirac-TDHF calculation in the GN model is the scalar potential S. Guided by the idea 
that the baryon potential has a kink-antikink shape as well as by covariance, we use as ansatz a rational function of 
an exponential U , 

A/^ _ ao + aiU + aiiU'^ 
^ V ^ 1 + biU + bnU^ ' ^ ' 

with 

C/ = exp{2?/7(a;-i)i)}, 7 ^ /J- — ^ - (5) 

V 1 — f 

y is a parameter governing the size and shape of the potential, v the baryon velocity. Note that for a static (w = 0) 
DHN baryon, the solution in ([2]) can be re-written in this form as 

S{x) = , U = e^y^ (6) 

and it is straightforward to boost this single baryon solution. Asymptotically, S must reach the value of the dynamical 
fcrmion mass, set to 1 by our choice of units, 

lim S = 1, lim S = 1. (7) 

This yields the conditions 

Ho = 1, On = &11. (8) 
By shifting the origin of the x or t axis we can rescale U and hence impose one further condition which we choose as 

an = 1, (9) 

leaving us with two unknown, real coefficients ai, &i. 
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Wc now turn to the Dirac-TDHF equation We use a chiral representation of the Dirac matrices, 7° = tTi,7^ = 
*o'2,75 — ^<^3, and introduce lightconc coordinates 

z — x — t, z — x + t, do — B — d, di—B + d. (11) 

U then becomes 

U = exp{y{Tj-^z + 7jz)} , (12) 

where 



r7 = e« = \/^ (13) 

turns out to be a more convenient variable than either the rapidity ^ or the velocity v. The Dirac equation assumes 
the simple form 

2idTjj2 = S4>i, 2idi'i ^ ~Sip2 (14) 

in terms of the chiral spinor components ipi = ipL,'4'2 = '4'R- 

In the next step, we try to solve Eqs. by an ansatz for the spinors. We assume from the outset that the 
potential is reflectionless and that the spinor has the same structure as S, multiplied by an exponential (plane wave) 
factor. In order to have a chance of solving the Dirac equation, the denominator will be taken to be identical to the 
one of S, whereas we admit arbitrary complex coefficients in the numerator. 

Consider the continuum spinors first. They are written as 

^k^ijd) p ^ cj = ±VT+^, (15) 



with 



Ni = cq + ciU + ciiU'^ , 

N2 = do + diU + dnU'^, (16) 



and V from ([T0| . The exponent can be rephrased in lightcone variables as 



kx - Lut = ^ ((z - ^ ] (17) 



with the spectral parameter 



(note the more symmetric definition of C as compared to Ref. jl7l|). Like z and z, gets simply rescaled under Lorentz 
transformations. The sign of uj is encoded in the sign of C, so that all results for continuum states hold for both signs 
of the energy. Notice also that it is not necessary to introduce two kinematical variables (lightcone momentum and 
lightcone energy), since only on-shell variables enter here. The asymptotics again puts constraints on the parameters 
of our ansatz for the spinors. Comparison with the solution of the free, massive (m ~ 1) Dirac equation. 



(0) - ^ f \ p»(C2-^/C)/ 



= -^^^ ( _\ ) eHC-^/C)/^ (19) 



yields 



Co = ■ .-, Cll=:JC0, 

do = 7^=, dn = Tdo. (20) 

VC' + l 
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Here, cq and dg are determined by the incident wave, cn and dn by the transmitted wave. T is the unitary transmission 
amplitude, \T\ = 1. The complex parameters ci,di,T have to be determined from the Dirac equation, together with 
the real parameters ai, bi in S. 

We now insert the ansatz for S and ipc into the Dirac equation. Since differentiation and multiplication yields again 
rational functions of the exponential U, we can equate powers of U, reducing the problem to an algebraic one. Owing 
to the use of lightcone variables, the calculation is straightforward, yielding the following results 

ai = {l-2y^)bu = w^^l-y\ 

w 

2{Z^-l + 2y') 



di 



-2(2ZV + 
iu{Z^ - 2iyZ - 1) ' 



Z'^ + 2iyZ - 1 

^ = Z^-2^Z-1 - (21) 

Z denotes the spectral parameter ^, boosted into the rest frame of the baryon. 

Let us consider the bound state spinors next. The transmission amplitude T has 2 poles in the complex Z-plane, 

Z^ -2iyZ -1 = {Z - Zi){Z + ZD, Zi=iy-w (22) 

corresponding to the positive and negative energy bound states at i? ±w in the rest frame of the baryon. Since the 
potential is transparent, the spinor may be regarded as a continuum spinor, analytically continued to complex ^. For 
positive energy spinors, this suggests the ansatz 



(1) \ 



with the exponent 



and 



= iy^{x~ vt) - wj{t - vx) = - [ Ciz ~ ] , Ci^V Zi (24) 



A/'i'^ = hU + hiU\ (25) 

Since 

le'^P = C/-\ (26) 

constant terms with eo,/o have to be omitted in the numerators (|25p for the sake of normalizability. If one inserts 
this ansatz into the Dirac equation, one can determine the unknown parameters algebraically except for the overall 
normalization with the result 

en = -ZiCi, 
/i = -Zirjei, 

hi - mi- (27) 
The parameter ei then follows from the normalization condition which we choose in the non-covariant form 

dx{ipHf^^ = 1, (28) 



yielding 



el = ^. (29) 

■qw 
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Before proceeding, we have to introduce one further element of hghtconc notation. One advantage of using C and rj 
rather than k and v is the fact that all square roots of the type Vl + A?, \/l — v'^ are eliminated. In order to avoid 
also w = y^l — y^, it is advantageous to combine y, w into the complex variable Zi introduced in (|22p. Using the fact 
that \Zi\ = 1, we get 

In this new notation the relevant single baryon quantities read 

2(Z4 + 1) 



01 = - 
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Z^{Zl-r\Y 

2(2^2^1^ - Z\ - 1) 
(Z2 + l)(Z-Zi)(ZZi + 1) 
2(^2^4 + Z2 - 2Zl) 



rfl = 

T = 



(Z2 + l)(Z-Zi)(ZZi + l) 
(Z+Zi)(ZZi -1) 
(Z-Zi)(ZZi + l)' 

^ (zl-v 
V \z! + i 



(31) 



All other coefficients are the same as in Eqs. ((20|) and ((27)) . A minor shortcoming of this notation is the fact that the 
unitarity of Zi {\Zi\ = 1) is not manifest. Thus for instance, the transmission amplitude T could also be represented 
as 

_ {Z-Z*){Z + Z,) 
{Z-Z,){Z + Zt) 

since l/Zi ~ Z*. This would show at once the poles of the positive {Z = Zi) and negative {Z ~ —Z*) energy bound 
states, as well as the fact that |r| = 1. To ease the notation of more complicated expressions later on however, we 
prefer not to use complex conjugate variables and will stick to the notation of Eqs. (|3ip . 

Notice that one could also generate the bound state spinor by analytic continuation from the continuum spinors 
(up to the normalization factor), rather than by ansatz. As one can easily check, it is sufficient to evaluate the residue 
of ip^ at the pole C = Ci ■ Finally, the negative energy bound state spinor can trivially be obtained from the positive 
energy one by the mapping "01 ~^ "01 j "02 — >■ —'02, so that there is no need to discuss it separately. 



B. Self-consistency and fermion density 

Having determined the free parameters in the ansatz for S and solved the Dirac equation, we proceed to verify 
self-consistency and derive the relation between the parameter y and the fermion number of the DHN baryon. To 
this end, we first need the scalar density ipTp for all occupied states. We decompose the scalar density for continuum 
states computed from the above spinors as follows. 



V'cV'C = (V'C'0c)i + (V'cV'c)2 

e + i 

2(^2(^2 + 1)2 



i^C^ch = ,.2 , .^;^2^2' .^^2 ^2^ (l-'g)^ (33) 



(C2 + l){Z^Zf - 1)(Z2 - Z2 

Positive and negative energy bound states on the other hand yield 



# = tUS^)(1-^)- (34) 



2 \Z( - 1 
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The scalar density of any occupied state can be written as a linear combination of two independent functions of (x, t), 
so that we are dealing with a type II solution of the TDHF problem. In the case of the kink, y = 1, Zi ^ i so that 
('0cV'c)2 a-nd the bound state scalar density vanish and the solution becomes type I. 
The self-consistency relation in the massless Gross-Neveu model is 

= -J^2S- (35) 

Rcnormalization of the coupling constant is performed with the help of the vacuum gap equation {S — 1), using a 
cutoff regularization 

dk dce + 1 ,,(0) 2c 

-A/2 27r 7i/A 27r 2C^ + I 



We find 



so that Eq. ([55]) becomes 



Ng 



^ ^ -InA (37) 



(V;V>=--lnA. (38) 

TT 

The first contribution to the chiral condensate from the continuum in Eq. p3p already gives self-consistency, 

^^(^C^c)i = --lnA. (39) 
i/A 27r 2C TT 



(V'cV'c)2 = -J— ^2 7 In , (40) 



Hence the 2nd contribution, 

/•OO 

io 27r 2(2 ^^^-^^^^^ 4^ V^?-l 

must cancel the contributions from the discrete states (occupation fractions h'±), 

2'Ki{v+ - v^) + = 0. (41) 

This determines the parameter y, given the occupation of the bound states. For the ground state baryons (y^ = 1) 
in particular, we recover the result of DHN, 



y = sin [-^ j . (42) 

In order to relate the occupation fractions v± of the bound states to the fermion number of the baryon, we still have 
to compute the fermion density, including the induced contribution from the Dirac sea. The density of the continuum 
states can be expressed through S as follows. 



/t/ {z^ + m^ + z^)zi ,^ 

^^^'^ ~ (C2 + 1)(Z2Z2 - 1)(Z2 - ZD ^ ^ 



(43) 



Subtracting the divergent density of the non-interacting Dirac sea, this yields the following contribution to the fermion 
density. 

Positive and negative energy bound states give the density 
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The total (subtracted) fermion density of the DHN baryon can therefore be expressed in terms of the bound state 
density as 

(^t^) = + -l)(^t^)(i\ (46) 

whereas the fermion number becomes 

Nf = N{iy+ + z/_ - 1). (47) 

If = 1 (ground state baryon), the induced fermion density from the continuum and the fermion density from the 
negative energy bound state cancel exactly, so that one can read off the full fermion density from the positive valence 
level alone. 

We finish this section with a comment on the massive GN model, i.e., Lagrangian ([1]) supplemented by a mass 
term — mo^/'V'. It is well known that the functional form of S of the DHN baryon also solves the massive GN model, 
the only difference being the relationship between y and i'± [H, [l^l . The self-consistency condition ((35)) has to be 
replaced by 

(^^) = -^(5-^o) = -^^+^ (48) 



with the "confinement parameter" 



The vacuum gap equation (|37p now reads 



so that becomes 



Trmo 



(49) 



-(7c + lnA) (50) 



(^^) =_^lnA+^(l-5) (51) 

TT TT 



It can be satisfied in the case of the DHN baryon because {ipijj) + (5/7r)lnA (which receives contributions from the 
continuum and the bound states) is proportional to (1 — 5*). 



III. BARYON-BARYON SCATTERING 



We now apply the ansatz method to a more difficult problem where the solution is not known — scattering of 
two DHN baryons with different baryon numbers (parameters 2/1,2/2)- To keep the number of parameters as small as 
possible, we shall work in the center-of-velocity frame where the baryon velocities are ±v. Since the calculation is 
fully covariant, we can transform the results into any other Lorentz frame afterwards. In analogy to the one baryon 
problem, we parameterize the scalar potential as a rational function of exponentials. The spinors are again taken to 
be exponentials times functions similar to S, where we always insist on keeping the same denominators. In the single 
baryon case, the asymptotic information which we used to reduce the number of parameters came from the vacuum. 
Now we will similarly exploit the asymptotic information from the incoming and outgoing baryons. This recursive 
way of proceeding greatly reduces the number of parameters which then have to be determined algebraically via the 
Dirac equation. 



A. Ansatz for scalar potential 

Our ansatz for the scalar potential is 



S-^ (52) 



with 



Af =1 + aiUi + a2U2 + anUl + ai2UiU2 + a22C/| + aii2C/f t/2 + ai22C/i?7| + aii22?7i^J7| 



(53) 
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V = l + biUi + b2U2 + biiUl + 512C/1C/2 + &22?/| + bii2UlU2 + &122C/1 C/| + bii22UlUl (54) 

and 

C/i = exp {2^17(0; - -yi)} , U2 ^ c^v {^y2l{x + vt)} (55) 

The Ui,U2 dependence is motivated by the product of the two baryon potentials, which we must recover when the 
scatterers are weU separated. In this sense, the ansatz is the minimal one having a chance of describing baryon-baryon 
scattering. Almost all of the 16 real parameters in S are determined by the asymptotic in- and out-statcs as follows. 
For t — >• —00, S reduces to 



l + aiUi +aiiUl 
c/2'A'o " l + hiUi+biiUl 



lim„g= \,Trr, ?rr\ (56) 



near incoming baryon 1 and to 



Ui^QO bii + bii2U2 + bii22U, 

near incoming baryon 2. For t — ^ oo, we get the limit 



c- an + aii2C^2 + aii22t/| 



,. c< 0-22 + ai22Ui + aii22U-^ 

hm S = ; — ; -py (58) 

C/2^oo ^22 + &I22C/1 + &1122C^i 



near outgoing baryon 1 and 



l + a2U2+ a22Ui 
uT-^o " 1 + 62 1/2 + 622 ui 



]^^S= \^,"[Z JZl (59) 



near outgoing baryon 2. These expressions should be matched to the single baryon scalar potential ([TU]). In this step, 
we have to account for the time delay occuring during the collision. To this end, we replace C/i — > Ui/ 612,1/2 -> 
U2/S21 (with real ^12,1^21) in the single baryon scalar potentials of the outgoing baryons. We only get a consistent 
parameterization for S12 — I/S21 ■— S and can now fix all coefficients in 5* except for 012, 612 and d as follows. 



ai 


= a{, 


02 = d. 


aii2 




ai22 = a{5, 


bi 


= bi 


62 = bi'5, 


6112 




6122 = b{S, 


ail 


~ aii22 


= 611 = &1122 = 


0-22 


= &22 = 


6\ 



(60) 

The superscripts on the one-baryon coefficients ai,6i on the right side of these equations refer to baryons / 

(parameters yi,v) and // (parameters y2, —v), respectively. As in the one-baryon case, it is advisable to use lightcone 
variables to keep the results for the spinors in a manageable form. We therefore write Ui,2 as 

Ui = exp{yi{r]^^z + r]z)} , 

U2 = cxp {y2ir]z + r]~^ z)} , (61) 
with 77 as defined in (|13p . and use the parameterization (i ~ 1,2) 

Zi = iyi - w^, \Zi\ = 1, 

The single baryon coefficients entering Eqs. ()60p then go over into 

2{Zt + 1) iZi 



a, = -- 



hi - 



zi{zi + iy 

u _ 2(Z| + 1) , _ 4Z2 

^2(^1 + 1)' ^ ^1 + 1- 

We will see later that b can actually be predicted on the basis of the single baryon input. This leaves us with only 
two unknown, real parameters ai2,&i2, to be determined from the Dirac-TDHF equation. 



^\ - v /'72 , °i - '72 , 1 ■ y^-^i 
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B. Ansatz for continuum spinors 

Anticipating that the mean field is refiectionless, we propose the following ansatz for the continuum spinor in TDHF, 

/ A/; \ p»(C5-^/c)/2 



Co 


- 

— Cq, 


Cl = 


do 


= 4, 


di = 


Co 


- C-, 


C2 = 


do 


= 4', 


d2 = 



C is the spectral parameter introduced in (|18p . 2? is the denominator of S, and the numerators are polynomials in 
Ui, U2 of the same degree as V, 

Ni=ca + ciUi + C2U2 + ciiUl + C12U1U2 + C22Ui + cn2U^U2 + ci22UiU^ + cn22U^Ul (65) 

AA2 = do + dif/i + d2U2 + diiUl + dr2UiU2 + d22t/| + dii2UlU2 + di22UiUi + dn22[/i [/|. (66) 

Most of the parameters are again determined by asymptotics. For the incoming baryon / and the outgoing baryon 
//, we can literally follow the treatment of S by letting C/1,2 — > 0, respectively, with the result 

Cl, Cii = Cii, 

d(, dii = d{^, 
ocj , C22 = d Cii, 

.8d{\ d22=<5M(- (67) 

We have taken into account the time delay for baryon // through the substitution IJ2 — > U2l^- For incoming baryon 
// and outgoing baryon /, we must account for the time delay of baryon / {U\ — SUi). In addition, here the incident 
continuum states have already been scattered by the other baryon. We therefore insert single baryon transmission 
amplitudes T^'^^ into the results analogous to S (for Ui,2 00), 

rrl II rp/ II rrl II 

cii = r Co , C112 = T Cl , C1122 = T Cj]^, 

dii - ^^d^^ dn2 = T'd[', du22 = T'dil 

C22 - T^Vc^, ci22=T"(5c(, cn22=T^Mi, 

d22 = T"6^d(,, di22^T"dd{, dii22 -T^dfi. (68) 

For convenience we list all single baryon parameters entering Eqs. (|67l68p . For baryon /, take over the results from 
Sec. Ill Al with the notation 

Z^Zi=i]C, Zi=iyi-wi, (69) 

i.e., 

r ^ 2{2Z]Z'i - Zf - 1) , 

' (Z2 + l)(Z,-Zi)(Z,Zi + l) °' 

d^ - ^ -T^d^ 

^ \/ C^ + 1 ' ^ 

, _ 2{Z]Zt + Zj 2Zl) 

" {Zf + l)iZi-Z,){ZjZ, + lf"- ^^"^ 
The analogous expressions for baryon // are 

Z ^ Zii = 7]~'^C, Z2 = iy2 - W2 (71) 



and 



// _ 'a II _ rpll II 

II 2(2Z|jZ| — Z| — 1) 



ji 



(Z| + 1)(Z„-Z2)(Z„Z2 + 1)'^° ' 

r]ii - ^ d"-T^^d^' 

VC + 1 

jII 2(Z|jZ| + Zfj - 2Z|) 

^1 ~ (Z| + l)(Zz,-Z2)(Z„Z2 + l)^° • ^'^^ 
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The transmission amplitudes can be inferred from Eqs. (j3ip . 



, ^ {Zi + Z,){ZiZ,-l) 
{Zi-Zi){ZiZi + l) 
^ {Zn + Z2){ZjjZ2-l) 



(Zii - Z2){ZnZ2 + l 

Some coefficients appear repeatedly in Eqs. ()67l68p . Due to the single baryon identities 



^-c^^ ci.^T^ci cii = T^^ci^ (74) 



and similar relations for c?'s. there is no conflict though. Eventually, only two complex parameters in the continuum 
spinor remain to be determined, ci2 and di2. Finally, we note that the transmission amplitude of the full continuum 
spinor factorizes, 

T = T^T'^. (75) 
This simply follows from the fact that we may evaluate it at a time where the baryons are well separated. 



C. Ansatz for bound state spinors 

It is sufficient to consider the two positive energy bound states. Along the lines discussed above, we set 



(i) 



where V is the denominator of S. Let us focus on the first bound state (asymptotically belonging to baryon /), since 
the 2nd one can simply be obtained by a relabeling of baryons / and //. As ansatz for the numerators, we choose, 

M''' - eiC/i + eiiC/f + ei2C/iC/2 + en2C/fC/2 + ei22C/iC/| + eii22C/fC/|, 

M'^ = /lC/i+/iiC/2+/i2C/lC/2 + /ll2C/fC/2 + /l22C/lf/| + /ll22C/fC/|. (77) 

The exponent is determined by the known bound state energy and kinematics, 

-FW = ^(Ci^--^), Ci=V-'Z^ (78) 

with Zl from Eq. dM]). Note that 

le'^'^'p^L/fi. (79) 



Normalizability then forces us to leave out all terms not containing Ui in ([77|. 

We now turn to the constraints from the asymptotic behavior of Since bound state 1 is attached to baryon / 
asymptotically, the only issues are the incoming and outgoing baryon /. For incoming baryon /, let U2 ^ and find 

ei = e{, en = e{i, /i - /ii = //i. (80) 

In the outgoing channel, let U2 —> 00. Here the bound state spinor 1 acquires a transmission amplitude from baryon 
// denoted as T/^. The asymptotic conditions therefore read 

ei22 = r/Me(, 61122 = T/^e(i, /122 = T/^'^//, fii22 ^ T^' f[,. (81) 

Let us collect the required one-baryon coefficients. 



e(i 


= -Ziei, 




// 


= -Zn]e{, 














1 


v\zU 


1 



(82) 
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How docs one compute a transmission amplitude for a bound state? Since the bound state spinor can be thought of 
as a continuum spinor with complex spectral parameter, we can find T/^ by analytic continuation. Take T^^ from 
Eq. ([73)) and replace the argument Zu by e~^^Zi = ri~'^Zi. e~^^ appears because we have to boost from the rest 
frame of baryon / (velocity v) to the rest frame of baryon // (velocity —v). In this way we arrive at 

1 (77-2Zi-Z2)(7?-2ZiZ2 + l)' ^ ' 

Since the spectral parameter is now complex, T/^ is not unitary. Its phase produces a phase shift of the bound 
state spinor, whereas its modulus gives rise to a time delay. This last observation gives us a clue for computing the 
time delay factor 6 introduced in the ansatz for S. Since the bound state spinor is moving along with the baryon 
asymptotically and the potential is quadratic in the spinor, S must satisfy the relation 

<5 = |r/^p. (84) 

Otherwise, there would be a spatial shift between the scalar potential and the bound state spinor of the scattered 
baryon. It is surprising that one can predict the time delay, a two-baryon scattering observable, on the basis of single 
baryon results alone. Evaluating ([M)) yields 

. ^ n ^ (ry^Zz + Zi)(r;2Zi + Z2)ir^^ZiZ2 - 1){7J^ - Z1Z2) 

d {v^Z2-Z,){tj^Z,-Z2){v^Z,Z2 + 1){v' + Z,Z2)- ^ ' 

Eqs. (|80H85p determine the parameters of the bound state spinor except for the 4 complex coefficients ei2, 6112, /12, /ii2- 

D. Nontrivial coefficients from the solution of the Dirac equation 

The ansatz for S and the spinors has been constrained as much as possible by single baryon results and asymptotics. 
This reduces the number of parameters from 16 to 2 real parameters in 5, from 18 to 2 complex parameters in the 
continuum spinor, and from 12 to 4 complex parameters in the bound state spinor. The remaining coefficients have 
to be determined by inserting everything into the Dirac equation (|T4)) and equating coeSicients of monomials C/f ?7™. 
The system is strongly overdetermined, so that the existence of a solution is non-trivial. We have indeed found a 
unique solution, confirming the correctness of our ansatz, with the following parameters: 

1. Scalar potential 

The non-trivial parameters 012, 612 are given by 

AZfZUZl + l){Zf + l)(7y« + 1) - 8{ZtZl + l){Zf + Z^Ty^ 



ai2 



ZiZ2iZl + 1){Z^ + l)d 
8Z1Z2 [2ZfZl{rj» + 1) - {Zf + \){Z\ + \)if\ 

{z-f + i){zl + i)d 



(86) 



with d from Eq. ((85 



2. Continuum spinor 

The parameters C12, di2 are found to be 

'^^•^12 , di2 
C12 = , , ai2 = . , (87) 

with the common denominator 

^12 - (Z2 + l)(Z| + l)(7?2ZiZ2 + l)(Zir/2-Z2)(r/2 + ZiZ2)(Zi-??2^2) 

(C - ^Z2){^Z2 + - z,){jiQZ^ + 1) 
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and the numerators 

ci2 = I^ZlZl [ZliZf + 1) + Z^iZl + l)?;i2] 

+ v')z',zi [(1 + zt){i + zl) + izfzic'] 

-CW {zl[,z\ + \){Zl + 1) + i^^Z\[Z\ + 1)(Z« + 1)] 

{[z\ + z\){z\z\ + 1) + C*^?^2'(^i + 1)(^2 + 1)] , 

di2 = -%Q^ZlZl{Zl{Zi + \)+r^''Zl[Z\ + \)\ 

+AZlZW{l + if) {C\Z\ + \){Z\ + 1) + 4Z?Z|] 

+4C2,74 [^2(^4 ^ i)(^8 ^ 1) ^ ^4^2(^8 ^ i)(^4 ^ 

-Sr/*^ {ZlZl{Z\ + 1)(Z| + 1) + C^(Zi4 + Zt){Z\Z\ + 1)] . (89) 



3. Bound state spinor 



The 4 relevant parameters are 

ei2 = 



2(ZiZ2 - 77^)(7?^Z2 + Z,)irj^Z!Zl + r^^Z^^ - 2Z|) ^, 
(Z| + l)rf 
2Zi(2Zi2z|-77''-r/'^Z|) j 
{Zl + l){Z,-r,^Z2){rf + Z,Z^f' 
2r;Zi(ZiZ2 - 7y2)(r;2Z2 + Z,){2fj^Z^Zl - 1 - Z\) j 



1 



(Z| + l)d 
2r,{ZlZl + - 2r?4z|) 

(Z| + l)(Zi-r;2Z2)(7y2 + ZiZ2)^ 

We reiterate that most of these expressions would look much more complicated if expressed in terms of conventional 
variables (fc, yi, j/2, w). Only the potential S and the time delay factor 5 are fairly simple in either representation. Since 
these are also the most interesting quantities, we have collected the results for S and 5 using conventional variables 
in the appendix for the convenience of the reader. 



E. Self-consistency and fermion density 

A crucial step in the TDHF approach still missing so far is the proof of self-consistency. We have to verify that the 
scalar condensate computed from our spinors reproduces the mean field S. We first evaluate the scalar density from 
the above spinors. The result for positive energy discrete states is 



1Z2 V 



with 

2Z, [2v^Zl{Zt + 1) - (778 + l)Zf{Zl + 1)] 



Xi 



{zl + \)d 



2Z, [2rfZf{Zl + (r;« + l)Zl{Zt + 1)] 
\2 = -■ (92) 

The negative energy states yield the opposite sign. In the case of the continuum states, we decompose the scalar 
density into 2 terms, following the same strategy as in the single baryon case, 

(V'cV'c)i = -^2-^*^' 

(^C^C)2 = J-<^'^^V^)(i)+^f (V^V)''^. (93) 
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f 1 2) 

The coefficients J-^ ' are (x, t)-independent. Since 3 independent functions of (x, t) are needed in order to represent 
the scalar condensate of an arbitrary single particle state, we are dealing here with a type III solution of the TDHF 
problem. Using the boosted spectral parameters Zj, Zjj introduced in ()69|71|) . the J"^^'^^ can be represented as 



i {e + i){zj-zi){zjzi-iy 
(2) _ ^iCzhizt - 1) 

[e + i)[zj,-zi)[zj,z'^-i) 



Integrating the first term in (|93p over dC, gives self-consistency, just like in the case of the single baryon. The integration 
over the 2nd term yields 



f 

Jo 



/o §^(V^C'/'c)2 = -^[(^^)(')lnZ^ + W)(^MnZ|]. (95) 

It cancels the contribution from the discrete states provided that the conditions 

27ri(i/,,+ - + lnZ.f 0, (j = l,2) (96) 

hold. These conditions are identical to the result for the single baryon, Eq. (|4T1) . confirming the self-consistency of 
the scattering solution. 

So far we have only dealt with the massless GN model. As discussed at the end of Sec. IIIBl the single baryon 
solution also solves the massive GN model, with a modified relationship between the parameter y and the occupation 
fractions h'±. As we have seen, a prerequisite for this to happen is the proportionality 

(?/;V) + -lnA-(l-S'). (97) 

In the baryon-baryon scattering case, this proportionality does not hold (no matter how one chooses the parameters), 
so that our ansatz does not lead to a solution of the massive GN model. Presumably, this refiects the fact that the 
massive GN model is not integrable, so that inelastic processes will contribute to baryon-baryon scattering if one 
switches on the bare mass (for a recent discussion of this issue, see Ref. [1^). 
The fermion density for the bound states can be computed from the spinors, 

^ (ij^ijf K (98) 

The analytical result is not very instructive and will not be given here (unlike in the single baryon case, it is not 
proportional to S*^ — 1). For the negative energy continuum states we find the subtracted density 

^JV^C-l = -^JV^-erp('\ (99) 
where the p^*^ are the densities of the discrete states. The C/^*'' are (x, <)-independent coefficients satisfying 

|4^ef = l, (^ = 1,2). (100) 

Hence the situation is again similar to the single baryon case. The total (subtracted) density can be computed from 
the bound state densities alone, 



Y.{>^^.+ +l^^,--l)p^'\ (101) 



and the total fermion number is just the sum of the fermion numbers of both baryons. For ground state baryons 
{vi^- = 1); the induced fermion density in the Dirac sea is cancelled exactly against the negative energy bound state 
contributions. Since the spinors have been given above and the analytical result for the density apparently cannot be 
simplified very much (unlike the scalar densities), we refrain from writing it down. 
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F. General Lorentz frame 



We have been working in the Lorentz frame where the baryons move with velocities ±v. Suppose we are interested 
in a general frame where the velocities are vi,V2- Since the Dirac-TDHF approach is covariant, we can transform the 
scalar potential and the spinors by a Lorentz boost. Consider the scalar potential S first. The velocity dependence 
enters in the exponentials Ui, U2, Eqs. (j6ip . in the time delay parameter S, Eq. (j85|) . and in the coefficients 012, 612, 
Eqs. (|86l) . We evidently have to replace Ui^2 by 



Ui = exp {yi{r]- ^ z + rjiz)} , rji 



1 



1 -Vi 

All coefficients entering S must be Lorentz scalars. Denoting the baryon 2-velocities by 



7» 



the only non-trivial Lorentz scalar available is 



U1U2 = 7172(1 - W1W2). 



In the frame where wi = — U2 = u, it reduces to 



(102) 



(103) 



(104) 



U1U2 



(105) 



We can then "covariantize" 5' by replacing 



2 U1U2 
V — > 



_ 7172(1 - viv-i) - 1 

U1U2 + 1 7172(1 - t'lW2) + 1 



where V12 is the (relativistic) relative velocity 



V12 



^12 



Vi - V2 



(106) 



(107) 



This rule is applicable to S or 5 written in conventional variables, as given in the appendix. As expected, the trans- 
formation to another Lorentz frame is more elegant in lightcone variables, where it reduces to the simple substitution 
rule 



V 



(108) 



This transforms the time delay factor 6 into the symmetric expression 

(772Z1 + r]iZ2){riiZi + 7/2 ^2) ('71^1^2 - '72)(??2^i^2 - m) 



S = 



{in Zi - rji Z2 ) (?/i Zi - ?72 ^2 ) {m Zi Z2 + ?72 ){ri2ZiZ2 + r}i)' 



(109) 



Likewise, the complexity of ai2, &12, Eqs. (|86p . docs not increase if one replaces ?/ by y/rjij r/2 . 

Finally, we remark that an even simpler form for 5 can be obtained by going back to Eq. (|83p for the bound state 
transmission amplitude. Let us transform T/^ into an arbitrary frame, using the substitution (|108p . 



//_ (cr + ci)(Ci*-C2) 
(cr-C2*)(Ci* + C2) 



(110) 



with Ci = 771 ^1 , C2 = 7?2 Z2 ■ This leads to the most compact expression for 5 in an arbitrary Lorentz frame which we 
could find. 



5 = 



(Cl+C2)(Cl-C2*) 



(Cl-C2)(Cl+C2*) 



(111) 
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The standard observable in Id scattering of solitons is the time delay, closely related to 6. Baryons / and // can be 
characterized by the following exponentials in the incoming and outgoing channels, 

C/i" = exp {2yi7i(x - v^t)} , U^' = S'^Ur, 

= exp {2y2j2{x ~ v^t)} , U^' = (112) 
The time delays Ati , Ai2 are introduced via 

Ur\x,t) = Ur{x,t~M,), (i = l,2) (113) 

and given by 

In (5 . In (5 , 
Ah = — , Ai2 = • 114 

2yi7iwi 2J/272W2 

So far, we have only discussed the transformation of S and 5 into a general Lorentz frame. The spinors can also 
easily be boosted in the lightcone approach, since the boost matrix e^^^/^ is diagonal. By contrast, the corresponding 
formulae in normal coordinates become exceedingly complicated. 



IV. LIMITING CASES AND ILLUSTRATIVE EXAMPLES 



We have presented above the general solution of baryon-baryon scattering in the large N limit of the GN model. 
In certain regions of parameter space, our solution can be compared to previous works. This gives us the opportunity 
to crosscheck our results against the existing literature. The special cases we arc aware of are the static limit 
(vi ~ V2 — 0), the non-relativistic limit (wi, ^'2, ?/2 ^ 1) and the kink limit (yi,i/2 — ^ 1), i-C-, the ultrarelativistic 
limit for the internal motion. We will also illustrate the full results with the help of some examples and make contact 
with the scattering of polarons and solitons in trans-polyacetylene in the present section. 



A. Static limit 



All static, transparent potentials of the Id Schrodinger equation have been constructed long ago by Kay and Moses 
Q. These Schrodinger potentials also enter in the construction of a complete set of static, transparent potentials of 
the Dirac equation, which are at the same time self-consistent potentials of the GN model 0, @ ■ Physically, they 
correspond to marginally bound n-baryon states whose mass is the sum of the constituent masses, independent of 
their separation. The n = 2 case can be related to the static limit [vi = W2 = 0) of baryon-baryon scattering. In the 
static case, one starts from the linear system of equations 

n 

^A,jVj = A,, (i==l,...,n) (115) 



where 



Ay=5y + ^^, A. =c.e-«-^ (116) 



The Dirac-HF potential is then given by [Tj 



s 



Comparing our results for u = with the static n = 2 solution, we get perfect agreement if we employ the following 
dictionary, 

yi ti, y2 = K2, 

Ui ^ A^e^y^'\ C/2 = A2e2^^^ 

cfwi \yi-V2 
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B. Non-relativistic limit 

Nogami and Warke [26| have constructed transparent potentials of the time-dependent, Id Schrodinger equation, 
thereby solving the non-relativistic TDHF problem for particles with (5-interactions. Equivalently, this amounts to 
finding solutions of the multicomponent nonlinear Schrodinger equation. We should be able to recover their results in 
the limit where both the internal motion of the fermions in the DHN baryons and their external motion as a whole are 
non-rclativistic, i.e., for yi,y2:Vi,V2 ^ 1. The non-relativistic, timc-dcpcndent construction of Ref. [26| generalizes 
the static one as follows: One starts from the linear system of equations 

±^+9. + u. = (119) 

with 

Ua = \/ Aa exp {kqiX + IK^^t] . (120) 

The Aa are real, the Kq complex parameters. The self-consistent potential can be written as 

V^2d^Y.<9o.- (121) 

a 

We have evaluated V for n = 2 and compared it with the non-relativistic limit of our S' in a frame where the baryon 
velocities are ui,U2. We treat yi,y2iVi,V2 as small quantities of order e and keep only the leading order term in 
(5 — 1), which is of order e^. Due to the fact that Nogami and Warke use units where m = 1/2 (to simplify the 
Schrodinger equation) whereas we use units where m = 1, we also have to rescale time t by a factor of 2 and compare 
their V with our 2(5 — 1). The two potentials then agree for the following choice of parameters, 

yi+ivi, K2= 2/2+^2, 
2yi, A2 = 2y2S, 

{yi + y2)^ + {vi - V2? ^^22) 

{yi - y2y + {vi - V2)^ ' 



Kl = 

Ai = 
5 = 



C. Kink limit 



In the limit yi — >■ l,i/2 1, the internal structure of the DHN baryon becomes ultrarclativistic. The baryons 
decouple into well separated kink and antikink. We should like to compare our result for S in this limit to the result 
for kink-antikink scattering [l^ [l^ , 

u cosh 27a; — cosh 27wt 
V cosh 27a; -I- cosh 2^vt 

Since the DHN baryon becomes infinitely extended in the kink limit, we have to shift our exponentials ?7i, U2 in such 
a way that only the scattering of the right edge of baryon / (an antikink) and the left edge of baryon // (a kink) 
survives in the limit j/1,2/2 ~^ 1- Working in the center-of- velocity frame, we first set 

yi = 2/2 = 1 - e. (124) 

The exponentials C/i, C/2 are replaced by 

Ui — A exp {27(2; — V i)}, C/2 = A-icxp {27(0; -Hwt)}. (125) 



For the choice 



one can then perform the limit e — > without encountering any singularity and S goes over into the negative of ()123p . 
The change of sign occurs because we get antikink-kink scattering, whereas (|123p holds for kink-antikink scattering. 
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FIG. 3: Scattering of a small [yi — 0.8) and a large (1 — j/2 ~ 10 ^) DHN baryon with velocities ±0.4. The self-consistent 
scalar potential S is shown for a range of {x,t) values in the vicinity of the collision. 
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FIG. 4: Same as Fig. |3l but total fermion density p shown. In the large baryon, the density has two peaks near the kink and 
antikink boundaries of the baryon. A negative time delay occurs during the collision. 



D. Examples 



In our first example, we consider scattering of a "small" baryon {yi = 0.8) and a "large" baryon (1 — y2 = 10"''), 
with velocities ±0.4. Fig. [3] shows the scalar potential as a function of x, for different time slices. The baryons cross 
each other with some (negative) time delay, as can also be seen from Fig. H] where the fermion densities are shown for 
the same scattering process. This should be contrasted to kink-antikink scattering where the scalar potentials first 
approach and then repel each other, the fermions being exchanged during the collision [isj . It is also interesting to 
watch the small baryon while it is crossing the large one in Fig. [31 In the initial and final state, the small baryon 
shows the ususal attractive potential well in the vacuum value S = 1, the large baryon exhibits well separated kink 
and antikinks with S = —1, i.e., the other degenerate vacuum, inbetween. During the collision, the small baryon 
moves as a seemingly repulsive potential bump on the S* = — 1 "floor" . This is just the chirally reflected DHN baryon 
in the other vacuum. To emphasize this aspect, we consider in our second example scattering of a small DHN baryon 
{yi = 0.8) on a CCGZ kink (1/2 = 1), with velocities ±0.2. In this case of baryon-kink scattering, the formula for S 
simplifies to the following expression (in conventional coordinates). 



S 



1 



l±^t7i±(5[/2 + [/? 



^V^J7iC/2 



with j/2 = 1 and 



(5 = 



l±i;^± 7/1(1-1;^) 
l±w^-yi(l-f2) 



(127) 



(128) 
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Comparing the small baryons in the first and last time slice of Fig. [SJ one sees clearly the sign flip of the potential 




FIG. 5: Scattering of a small DHN baryon {y = 0.8) on a CCGZ kink (y — 1), at velocities ±0.2. Notice the chiral reflection 
which the small baryon undergoes during the collision. 

characteristic for the discrete 75 transformation during the collision. 

The last few examples are applications of the GN model to condensed matter physics. As is well known, the 
Su-Schrieffer-Heeger theory (23] of trans-polyacetylene admits a continuum approximation [28j closely related to the 
massless GN model. At the mean field level, the correspondence with the large N GN model is perfect, although N 
in the condensed matter case is only 2, the number of electron spin components. DHN baryons become polarons, 
CCGZ kinks become solitons in conducting polymer language [29| . The polaron has one electron in the upper bound 
state, correponding to the DHN parameter y ~ sin7r/4. For polaron-polaron scattering with relative velocity v, we 
thus have to choose yi — 1/2 = 1/a/2 and find 6 — 



S 



1 + 2V2U1 



(129) 



Polaron-soliton scattering on the other hand corresponds to the parameters ±v and yi 
simplifies to 

^ 1~SU2 + Ul - UIU2 



1 + 2V2f/i + 5U2 + Ul + 2V2VSU1U2 + UIU2 



1/V2,y2 = 1, where S 
(130) 



with 



' {I + v^)^/2 + {1 ~ v^) 
(l + w2)V2-(l-w2) 



(131) 



Finally, soliton-soliton scattering can be inferred from Eq. (|123p in Sec. IIII CI 



V. SUMMARY AND OUTLOOK 



In this paper we have presented a complete description of the scattering of two DHN baryons, each with its own 
internal structure of localized mean fields and bound fermions, in the large TV limit of the Gross-Neveu (GN) model. 
We have used the time-dependent Hartree-Fock (TDHF) method, thereby realizing an explicit example of the picture 
envisioned by Witten in his seminal paper on baryons in the 1/N expansion [l3| . Static baryons of this model have 
been known since the early work of DHN, but progress on time-dependent, self-consistent solutions has been slow until 
very recently. The only cases which had been understood so far are the non-relativistic limit, where the model reduces 
to the multicomponent nonlinear Schrodinger equation, and the ultrarelativistic limit of kinks and antikinks. Here 
we have presented the most general solution, valid for 2 DHN baryons with arbitrary fermion numbers and velocities. 
This contains in particular scattering of polarons and solitons in the continuum limit of conducting polymers like 
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trans-polyacctylene. The solution displays characteristic features expected from an integrable model: Purely elastic 
scattering, vanishing reflection coefficients, baryons crossing each other with a time delay. The complete time evolution, 
including the non- asymptotic region, is available in closed analytical form through the TDHF spinors. The calculation 
is significantly more involved than in the kink-antikink case, but remains manageable owing to the systematic use of 
lightcone coordinates and variables with their simple transformation properties under Lorentz boosts. 

More interesting than the specific results for baryon-baryon scattering is perhaps the method used to solve this 
problem. Since the scattering solution is of type III, the TDHF problem is mathematically equivalent to solving 
3 coupled, nonlinear partial differential equations. It does not seem very likely that one can guess the analytical 
solution of such a complicated problem. Nevertheless, the method of ansatz has proven very effective. The basic 
idea was to use a simultaneous ansatz for the TDHF potential and the single particle spinors. This ansatz is based 
on very few heuristic rules: The scalar potential is assumed to be a rational function of 2 exponentials, one for 
each baryon. The degree of the polynomials in the numerator and denominator is taken from the product of two 
independent single baryon potentials. The coefficients are determined to a large extent by the known baryons in the 
incoming and outgoing states. The spinors are also assumed to be rational functions of exponentials times a plane 
wave factor (for continuum states) or its analytic continuation to complex spectral parameter (for bound states). The 
denominators are taken to be the same in S and all the spinors, the numerators have different coefficients fixed largely 
by asymptotics. When the dust has settled down, only a handful of coefficients remain to be determined by solving the 
Dirac equation, which is now turned into a simple algebraic problem. Having found a solution of the Dirac equation, 
we still have to check self-consistency of the whole approach. This was successful for the massless GN model but not 
for the massive one, although single DHN baryons solve both of these models. 

Our method can be generalized in a natural way to other, even more complicated, baryon scattering problems. We 
have already used it to reproduce the DHN breather and to solve a dynamical 3-baryon problem. This will be reported 
elsewhere. One could also consider the scattering of a DHN baryon from a lump of finite-density matter, represented 
as an array of kinks and antikinks, or even the scattering of two such lumps. The question then arises whether one 
can find the general mu/ti-baryon solution of the TDHF equation in closed form, similar to what has already been 
achieved in the static limit, in the non-relativistic limit or in the kink limit of the GN model. The two-baryon solution 
discussed here has many encouraging features which suggest that such a generalization might actually exist. If one 
can find it, it will be interesting to see whether it coincides with some other solved problem in nonlinear mathematical 
physics, like in the case of the sinh- Gordon theory, or whether it leads to a new class of exactly solvable problems. 
We hope that having worked out the 2-baryon problem in full detail, and having cast the solution into the simplest 
form, will help to answer this question in the future. 
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APPENDIX: SCALAR POTENTIAL IN CONVENTIONAL VARIABLES 

The self-consistent scalar potential for baryon-baryon scattering has the general form introduced in Eqs. 
Here we list all the parameters in conventional variables, complementing the lightcone variables used in the main text. 
Parameters determined asymptotically by onc-baryon data, 

2(1 -2^;^) 2(1 - 2yl) ^ 
ai = , 02 = 0, 

2(1 - 2yi) _ 2(l-2y?) 
aii2 — , ai22 — 0, 

2 2 
bi = — , b2 = — S, 

Wi ' W2 

h 2 2 

Oll2 — , 0i22 — 0, 

W2 Wi 

an — aii22 = &ii = ^1122 — 1, 

a22 = 622 = <5', (132) 
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with Wi — \/i~-yf and 



n ^ 4v^ + il- v^fivl + yj) + 2(1 - v^)yiy2 

d Av^ + (1 - t;2)2(y2 + y2) _ 2(1 - ^;4)yiy2 " ^ ''''^ 



Parameters determined via the Dirac equation, 

4[4i.2 + (1 - 10i>2 + i;4)(y2 + y2) _ 2(1 - v^)^{yf + + 2(1 + 6v^ + v^)yly^] 



ai2 



WiW2d 

A[4v^ + il-v^riyl+yl-2ylyl)] 



h2 = ' ' ' (134) 

Wi'W2d 

with d from Eq. 

In a Lorentz frame where the baryons move with velocities wi,U2, replace Ui,U2 in Eqs. (|55p by 



[/j = exp{22/i7i(x - Wit)) , 7^ = r— = t (* = 1:2), (135) 

V 1 - 



and replace w by the relative velocity vi2, Eq. (|107p . in S, ai2, 612- 
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